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an  ionic -dipol ar  attraction.  In  this  paper  we  consider  composite 
functions  which  consist  variously  of  an  exponential  form 
of  repulsion  or  an  inverse  R-12  repulsion,  and  attractions 
which  are  expressed  in  terms  of  the  basic  ionic -dipolar  terms 
together  with  additional  terms  which  depend  upon  the  polariz¬ 
abilities  of  the  ion  and  the  solvent.  No  more  than  two 
adjustable  parameters  are  used.  The  force  constants  for  the 
far  infrared-active  vibrations  of  the  caged  ions  are  known. 

With  the  use  of  the  equilibrium  condition  for  the  cage  of  solvent 
the  values  of  the  parameters  can  be  determined*,  Although  all 
of  the  functions  examined  show  the  effect  of  tpe  polarizability 
of  the  ion  on  the  force  constants,  a  potential'  v^ich  consists 
of  independent  R'6  and  R'l^  terms  together  with  other 
experimentally  fixed  ionic-dipolar  terms  shows  most  clearly  the 
effect  of  polarizability.  We  also  find,  however,  that  the 
consideration  of  ionic  polarizability  alone  is  inadequate  to 
explain  in  a  simple  manner  the  variation  of  the  force  constants 
for  the  alkali  metal  cations. 
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Abstract 


The  purpose  of  this  paper  is  to  explore  the  use  of  several 
model  potential  energy  functions  in  order  to  try  better  to  under¬ 
stand  the  nature  of  the  forces  which  operate  between  simple  monatomic 
cations  and  the  molecules  of  solvent  which  surround  them.  This  work 
builds  upon  and  extends  a  recently  reported  initial  treatment.  In 
that  report,  the  potential  energy  was  modelled  as  the  sum  of  an 
exponential  (Born-type)  repulsion  and  an  ionic-dipolar  attraction. 

In  this  paper  we  consider  composite  functions  which  consist  variously 
of  an  exponential  form  of  repulsion  or  an  inverse  R-12  repulsion, 
and  attractions  which  are  expressed  in  terms  of  the  basic  ionic- 
dipolar  terms  together  with  additional  terms  which  depend  upon  the 
polarizabilities  of  the  ion  and  the  solvent.  No  more  than  two 
adjustable  parameters  are  used.  The  force  constants  for  the  far 
infrared-active  vibrations  of  the  caged  ions  are  known.  With  the  use 
of  the  equilibrium  condition  for  the  cage  of  solvent,  the  values  of 
the  parameters  can  be  determined.  Although  all  of  the  functions 

examined  show  the  effect  of  the  polarizability  of  the  ion  on  the  force 

~6  ~12 

constants,  a  potential  which  consists  of  independent  R  and  R 
terms  together  with  other  experimentally  fixed  ionic -dipolar  terms 
shows  most  clearly  the  effect  of  polarizability.  We  also  find, 
however,  that  the  consideration  of  ionic  polarizability  alone  is 
inadequate  to  explain  in  a  simple  manner  the  variation  of  the  force 
constants  for  the  alkali  metal  cations. 
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1 .  Introduct ion 

It  is  generally  assumed  that  for  most  cations  in  solution, 
and  even  for  those  for  which  no  effective  covalent lv  bonded  associa¬ 
tions  form,  a  discrete  and  persistent  structure  of  solvent  molecules 
accompany  a  centrally  located  ion.  One  can  propose  significant, 
discrete  structures  which  may  contribute  to  chemical  reactivity  or 
transport  processes  which  take  place  in  solution.  Thus,  for  a 
variety  of  reasons  there  is  a  need  to  try  to  understand  the  nature 
of  the  forces  which  hind  solvent  to  ion  in  such  an  impermanent 
environment  as  that  of  a  solution. 

Tt  is  the  purpose  of  this  paper  to  present  several  forms  of 
potential  energy  function  for  the  simple  solvated  alkali  metal 
cation.  An  ultimate  objective  of  this  kind  of  investigation 
should  he  to  specify  an  accurate  form  of  the  function.  No  claim  is 
even  hinted  that  we  are  anvwhere  near  that  goal.  On  the  contrary, 
here  only  short  -  ranged ,  intermed late- -but  necessary- -oh j ect ives  arc 
set.  Specifically,  we  fit  parameters  in  these  model  potential 
functions  with  available  spectroscopic  data.  At  this  stage  it  is 
necessary  to  see  if  it  is  possible  to  choose  one  form  rather  than  another 
as  being  more  "accurate"  with  respect  to  physical  interpretation. 

As  we  amply  illustrate  in  this  paper,  the  stage  has  not  been  reached 
where  it  is  possible  to  assay  the  numerical  accuracy  of  any  particular 
form  of  potential  energy  function.  For  that  assay,  different 
experiments  with  new  sets  of  data  are  required. 

As  we  indicated  in  a  recent  paper, ^  the  solvated  alkali  metal 
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cations  are  ideal  subjects  to  study  because  of  the  relatively 

uncomplicated  nature  of  the  electronic  structure  of  the  centrally 

located  ions.  Thus,  one  suspects,  the  predominant  forces  of  inter' 

action  are  elect rostat ic  with  strong  repulsions  entering  only 

on  very  close  approach  between  the  ion  and  its  solvent.  The 

process  of  modelling  the  appropriate  potential  energy  function  is 

thereby  simplified.  Experimental  evidence  in  the  form  of  the 

2  .  - 

far  infrared-active  vibrations  of  the  alkali  cations  definite]}' 
suggests  that  the  ions  occupy  some  form  of  cage  of  solvent  and  that 
the  cage  endures  for  a  sufficient  length  of  time  for  the  vibrations 
to  be  observed. 

In  the  previous  paper, ^  we  modelled  the  interactions  between 
the  centrally  located  ion  in  a  cage  of  solvent  and  the  molecules 
of  solvent  which  make  up  the  cage  as  the  simple  sum  of  an  ion-dipole 
interaction  and  an  exponential,  Born-type  repulsion.  An  equilibrium 
condition  for  the  single  ion-solvent  couple  was  used  to  fix  values 
for  the  adjustable  parameters  in  the  Born  repulsion.  Most  importantl 
we  showed  that  because  the  ion  oscillates  about  a  position  of 
equilibrium  which  also  is  the  centre  of  symmetry,  the  force  constant 
does  not  depend  directly  upon  any  electrostatic  ionic -dipolar  terms. 
The  contribution  of  the  electrostatic  interaction  between  the  ion  and 
a  regular,  crystalline  cage  of  dipoles  vanishes  for  the  force 
constant  by  reasons  of  symmetry.  The  form  of  the  force  constant 
which  we  obtained  indicated  a  direct  dependence  on  the  forces  of 
repulsion  which  operate. 

We  also  showed  in  the  previous  paper^  that  the  force  constants 
for  the  vibrations  of  the  alkali  cations  in  their  cages  of  solvent 
range  from  a  high  value  for  lithium  to  a  low  value  for  potassium. 
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liowever,  from  potassium  to  caesium,  the  force  constants  again  increase. 
We  indicated  that  this  effect  could  be  attributed  to  the  increasing 
polarizability  of  the  electrons  in  the  ions  as  the  group  is 
transcended  from  lithium  to  caesium.  Our  arguments  in  that  pre¬ 
liminary  report, 1  however,  were  only  qualitative.  An  aim  of  this 
paper  is  to  attempt  to  make  the  dependence  of  the  force  constant 
on  the  polarizability  of  an  ion  explicit  and  quantitative.  We 
want  to  see  if  the  trend  in  the  force  constants  for  the  group  of 
the  alkali  metal  cations  follows  the  trend  in  the  increase  in  the 
observed  polarizabilities. 

We  consider  the  following  potential  energy  functions.  First, 
we  re-examine  the  composite  function  which  consists  in  part  of  the 
Born  exponential  type  of  repulsion.  Several  years  ago,  Woodcock^ 
considered  a  generalized  form  of  exponent ial - type  of  repulsion 
in  connection  with  a  study  of  the  alkali  halide  solids.  A  Woodcock- 
form  of  repulsion  is  used  in  the  second  example  function.  An  advantage 
of  the  Woodcock  generalized  repulsion  is  the  fact  that  the  "hardness" 
of  the  repulsion  can  be  adjusted  easily.  Finally,  we  consider 
composite  functions  which  are  constructed  of  Lennard- Jones  and  Mie 
(or  decoupled  Fennard- Jones)  potentials.  We  will  demonstrate 
that  the  use  of  the  Lennard- Jones  potential  in  a  Stockmayer- f orm 
of  potential  leads  to  pathological  conditions.  In  particular, 
one  is  not  free  to  choose  values  of  the  radii  of  solvation  for  the 
ions  at  will.  The  use  of  the  Mie  potential,  on  the  other  hand, 
shows  a  remarkable  R  ^  repulsion  for  the  interaction  of  the  ion 
with  the  solvent  both  for  lithium  and  sodium. 

Before  we  present  the  analyses  of  the  various  potential  forms 
in  detail,  it  is  necessary  to  consider  some  procedura 1  and  mathematical 


preliminary  operations.  These  operations  are  outlined  in  the  next 
sect  ion . 

2.  Method  of  approach 

In  this  paper  we  consider  the  same  systems  as  we  considered 
previously,1  viz.  ,  the  solvated  alkali  metal  cations.  As  indicated 
above,  here  we  consider  more  complicated,  and  one  presumes,  more 
realistic,  potential  energy  functions.  These  potential  energy 
functions  are  constructed  as  the  sums  of  terms:  repulsions  and 
attractions.  In  particular,  we  consider  various  ionic-dipolar 
attractions,  dipolar  repulsions  between  molecules  of  solvent  in 
the  first  shell  of  solvation,  dipolar- induced  dipolar  attractions 
between  solvent  molecules  and  between  the  solvent  and  the  ion,  and 
finally  in  some  cases,  induced  dipole-dipolar  interactions  between 
the  ion  and  the  solvent.  In  all  cases,  the  polarizability  of  the 
solvent,  as  found  with  the  use  of  the  Lorenz - Lorent z^  equation,  is 
used.  In  some  of  the  cases,  the  polarizability  of  the  ion  is  used; 
the  values  arc  those  given  in  Smyth's  book.^  The  parameters  which  are 
adjustable  and  require  fitting  to  ex  perimcntal  data  are  contained 
in  the  exponential  repulsions  or  the  Lcnnard- Jones  (or  Mic)  components. 
Throughout,  we  allow  for  only  two  adjustable  parameters.  These 
two  parameters  can  he  fitted  to  values  of  the  force  constants  for 
the  far  infrared-active  vibrations  of  the  ion.  The  expression  for 
the  force  constant  and  the  condition  for  equilibrium  yield  two 
expressions  to  use  to  solve  for  the  two  parameters. 

l;or  the  solvated  ionic  systems,  there  are  only  limited  un- 
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ambiguous  experimental  data  available  which  can  be  used  to  determine 
the  parameters  in  any  model  potential  energy  functions.  Indeed, 
as  long  as  the  question  of  the  coordination  number  of  the  solvent 
about  any  one  of  the  alkali  metal  cations  remains  open,  there  are 
no  data  which  are  certain.  We  assume  in  this  paper  that  the 
coordination  number  of  solvent  for  lithium,  sodium,  and  potassium 
is  four.  The  coordination  number  of  solvent  for  rubidium  and  caesium 
is  most  likely  to  be  six.  It  is  then  possible  to  establish  the 
values  for  the  force  constants  for  the  far  infrared-active  vibrations 
of  these  ions.'* 

Data  which  are  used  in  the  following  section  to  establish 
values  of  the  adjustable  parameters  within  each  model  potential 
energy  function  are  given  in  Table  1.  The  values  of  the  radii 
of  solvation  are  those  used  previously'*  and  are  found  with  the 
use  of  the  Stearn-Eyring^  model  as  applied  by  Abraham  and  Liszi. 

The  polarizability  of  the  solvent  (taken  to  be  dimethyl  sulphoxidc, 
DMSO)  is  found  with  the  use  of  the  horenz-Lorentz^  equation 

2—1  V  =  4ira/3  (2.1) 

n2  +  2  m 

in  which  n  is  the  index  of  refraction,  V  is  the  molar  volume,  and 
a  is  the  polarizability. 
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a.  Formulae  for  the  force  constant 

As  we  showed  in  the  previous  paper, ^  when  considering  the  far 
i n t rared -act i ve  vibrations  of  a  solvated  ion,  it  is  essential  to 
develop  expressions  for  the  force  constants  in  an  appropriate  form. 

The  vibrations  are  known  experimentally-'  to  involve  the  oscillation 
of  the  ion  about  a  position  of  equilibrium  at  the  centre  of  a 
cage  of  solvent.  IVe  assume  the  cage  is  regular  and  crystalline. 

Thus,  the  position  of  equilibrium  is  also  the  centre  of  symmetry. 

It  is  necessary,  therefore,  to  obtain  the  force  constant  as  the 
second  order  coefficient  in  a  Taylor  series  expansion  about  the 
centre  of  symmetry. 

For  the  electrostatic  interaction  between  the  ion  and  its 
solvent  dipoles,  the  Carl son-Rushbrooke®  expansion  can  be  used.  Tor 
cages  of  solvent  which  are  distributed  about  the  ion  with  cubic 
symmetry,  these  electrostatic  terms  do  not  contribute  to  the  magnitude 
of  the  force  constant.  For  other  potential  functions  an  integral  form 
of  the  Taylor  scries  has  proved  useful.  This  form  was  presented 
to  second  order  in  the  previous  nancr.^  Recently,  we  have  generalized 
the  result  to  any  order.  ^,10  The  essential  formulae  are  listed  below. 

Frequently,  one  encounters  potential  energy  functions  which 
do  not  easily  yield  Fourier  transforms.  The  need  for  the  Fourier 
transform  of  a  function  is  essential  in  evaluating  the  integral 
form  of  the  force  constant - -see  below.  When  the  Fourier  transform 
cannot  he  found,  vet  the  function  is  otherwise  well-behaved, 
it  is  neccssarv  to  revert  to  a  differential  form.  Therefore,  we 
also  present  below  a  differential  form  which  is  useful  in  evaluating 
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tile  force  constant  associated  with  the  Wookcock  generalised  exponen¬ 
tial  form. 

The  force  constant  for  the  vibration  of  an  ion  about  its 
equilibrium  position  at  the  centre  of  a  cage  of  solvent  is  obtained 

as  follows.  The  formula  for  the  second  order  contribution  to  the 

'r  i  •  -  1,9 

fa v lor  series  is 


fa 1 


vr)2c.(R)  = 
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. . c  +  2 , 
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(2.2) 


where 


i2£(R) 


l 


,00 


( 2 if)  3 J  o 


elk 


k4f(k).if  (kin 


(2.3) 


and  j  (x)  is  the  spherical  Resscl  function  of  the  first  kind  S'* 

a  a 

In  eqn  (2.2)  P  (x)  is  the  Legendre  polynomial . ^  The  argument  r*R 
of  the  Legendre  polynomial  is  the  scalar  product  of  the 

A  A 

unit  vectors  r  and  R.  The  vector  R  is  directed  from  the  molecule 
of  solvent  to  the  centre  of  symmetry  and  the  vector  r  expresses 
an  excursion  away  from  that  point.  The  function  f(k)  is  defined 


fOO 


(4  T  ) 


3/2 


dr  r2C(r) j  0  Or) 


(2.4) 


when  G(r)  is  a  scalar  function.  The  coefficient  lias  two  values 

1/3  and  2/3  according  as  l  =  0  or  2. 

For  a  system  of  solvation  which  is  regular  and  crystalline 
in  either  a  tetrahedral  or  octahedral  sense,  it  is  easily  shown 
that  the  the  9-2  terms  vanish  when  summed  over  the  sources. 
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Thus ,  the  contribution  to  the  force  constant  from  an  individual 
term  in  the  potential,  (Hr),  is 


ck  =  -  c/4T 1 2  o  (Re) 


•  OO 

dk  k4f(k)jo(kR0) 

0 


(2.5) 


where  c  is  the  coord inat i on  number  for  the  cage  of  solvent. 
If  the  function  G (r)  is  the  exponential  Born  repulsion 


exp 


B  exp [ - (R- R o ) / p] 


then  this  term  will  contribute 


(2.61 


\xn  =  —  (1  -  2p/R0)  (2.71 

exp  3p2 

to  the  total  force  constant.  ^  On  the  other  hand,  for  a  term  in 
the  potential  of  the  form  R  m,  a  contribution  of  the  form 


c 


k 


m 


m (m- 1) 
3R0m+2 


(2.8) 


is  found. ^ 

liquation  (2.4)  defines  the  radial  component  of  the  Fourier 
transform  of  a  function  G(r).  As  indicated,  frequently  one  encounters 
well-behaved  functions  which  nevertheless  fail  to  yield  easily 
analytical  expressions  for  f(k).  In  some  instances,  although  it 
is  possible  to  find  an  analytical  form  for  f(k),  it  is  not  possible 
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easily  to  find  the  integral  (2,3).  On  the  other  hand,  such  functions 
may  yield  readily  to  differential  operators.  One  function  which 

satisfies  tins  classification  is  the  Woodcock  generalized  exponential 

,  .  3  a 

repuls l on 

*W(R)  =  BR~n'cxp  [  -  (Rn-Ro*)/ pn]  •  (2.9) 


We  now  develop  the  differential  form  for  the  force  constant  term 

(2.5).  The  complete  expression  for  the  second  order  term  in  the 

1  9 

favlor  series  is  given  as  ’ 


t  2  ( r  ,  R ) 


jl 2  o (R)  +  jP2 (r-R)I22  (R) 


(2.10) 


for  a  function 


G(r)  =  /47  Yo0(r)F'(r) 


(2.11) 

.a 

in  which  Y0o(r)  is  the  zeroth-order  spherical  harmonic  function. 

The  problem  is  to  convert  l2o(R)  and  l22(R)  to  differential  forms. 

First,  for  l2o(R)  v?e  write  from  the  definition  of  the  integral, 


eqn  (2.3) 
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Into  the  integral  on  the  second  line  of  eqn  (2.1?)  wc  substitute 
the  definition  of  f(k).  Changing  the  order  of  integration  gives 
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f  dr  r^2F(r)6(r-R) 
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=  4tt5/,2F(R) 


(2.13) 


where  <5(r)  is  the  Dirac  delta  function.  Substitution  of  this  result 
into  the  equation  for  l2o(R)  gives 


Iao(R) 


1 

/4rf  R 


3 2  F 
3R2 


(2.14) 


For  the  evaluation  of  the  second  term  I22(R),  we  use  the 
lowering  operator  for  the  spherical  Ressel  function  twice.  * 


jn+l(x)  =  [n/x  -  d/dx] j n (x) .  (2.15) 

Thus,  by  means  of  an  operation  similar  to  the  above  to  get  (2.14),  wc 
find 


1 2  2 (R)  = 


/4it 
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R3R 


F  (R) 


(2.1b) 


The  collected  result  is 
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The  form  of  (2.17)  is  the  same  as  (2.10)  as  far  as  any  angular 
dependencies  are  concerned.  Therefore,  for  collections  of  sources 
which  are  distributed  with  cubic  symmetry,  tire  l\(r*R)  terms  vanish 
when  summed  over  all  the  sources.  The  application  of  the  result 
to  the  Woodcock  generalized  exponential  potential,  eqn  (2.9), 
yields 


cr 


W  l  in ,  n  ) 


r r2 


m(m-i)  -  n (n- 2m+l ) (R0/p)  +  n2(R0/p) 


cn 


(2.18) 


The  contribution  to  the  force  constant  for  a  Born  exponential 
repulsion,  eqn  (2.7),  follows  immediately  from  (2.18)  when  m  =  0 
and  n  =  1 . 


b.  Li  Lect  ros  tat  ic  interactions 


The  interactions  between  the  ion  and  the  solvent  dipoles  and 
between  the  dipoles  of  the  solvent  themselves  will  have  a  constant 
representation  in  each  of  the  potential  energy  functions.  The 
inclusion  or  omission  of  induced  dipolar  interactions  for  the  ion 
varies  with  the  type  of  repulsion  considered;  the  Lennard-Joncs  or 
Mie  potentials,  for  example,  implicitly  include  the  induced  dipolar 
interaction.  Here,  we  consider  a  general  form  to  be  used  through¬ 
out  the  disucssion  in  the  next  section. 
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From  Hottchcr  we  can  write  the  potential  energy  for  the 
interaction  between  two  identical  dipoles  of  solvent  as 


dd 


A 


b  1  #  b  2 


3A2  (Vji's)(W2-s) 

A7A  4  5 
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A  1  s  3 


3  cr2)2' 

A  l  A  4  s  5 


(2.10) 


in  which 


A  ■  1  -  n^^-(3/8)3 
n  ,,  6 


(2.20) 


and  p  is  the  magnitude  of  the  dipole  moment.  The  polarizability 
of  the  solvent  is  as •  Finally,  s  is  the  distance  between  the  centres 
of  gravity  of  the  dipoles. 

The  contribution  from  the  interaction  between  the  ion  and 
a  dipole  is 


(j) 


id 


1_ 
A  4 


(b  -R) 
Ze  - 
R3 


A  i  A  t 


(us-R) 
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Ze 2  a 
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2A4R4 


(2.21) 


in  which  R  is  the  distance  between  the  ion  and  the  centre  of  gravity 
of  the  dipole.  The  polarizability  of  the  ion  is  a..  Now,  in  A  , 
one  of  the  polarizabilities  is  that  of  the  ion  and  the  other  is  that 
of  the  solvent,  as. 

It  is  not  difficult  to  show  that  for  the  distances  applicable 
to  the  solvated  ionic  species,  the  A  all  are  essentially  unitv. 

In  the  following,  we  will  use  the  simpler  forms  in  which  all  A  =  1. 

The  distance  R  between  the  ion  and  the  dipole  is  the  parameter 


of  interest.  Therefore,  it  is  necessary  to  express  the  distance 
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s  between  the  molecules  of  solvent  in  terms  of  R.  for  the  regular 
tetrahedron  and  octahedron,  this  is  a  simple  process.  Furthermore, 
it  is  possible  to  relate  the  orientations  of  the  dipoles  to  the 
unit  vector  R.  for  the  \f'n  molecule.  If  we  assume  an  orientation 
which  yields  the  optimum  interaction  between  the  ion  and  the  solvent 
dipole,  then  the  dipole  vectors  can  be  written  as 
iii  =  uRi 
U  2  =  M  R  2 

for  a  pair  of  dipoles,  and 
s  =  I  R  a  -  R  2  |  • 

For  any  pair  of  dipoles, 
y  i  •  n  2  =  m  2  c  o  s  0 

=  u2Ri -R2 , 

and 

U  i  •  s  =  mRi  •  (P2-R1)  =  ViR(cosO-l) 
u 2  *s  =  -  uR(cos9-l )  . 

In  calculating  the  total  contribution  to  the  potential  energy 
function  due  to  the  interactions  between  all  the  dipoles  in  a 
tetrahedron,  we  multiply  the  basic  result  for  a  single  pair  by  b 
(i.e.,  n(n-l)/2  for  n  identical  species).  On  the  other  hand,  for 
the  octahedron,  we  multiply  the  basic  result  by  12  (i.e.,  n(n-l)/2 
-3  which  eliminates  the  three  pairs  which  are  separated  by  the 
centrally  located  ion). 

The  specific  forms  for  the  d i pole -d i po 1 e  and  ion-dipole 
energies  arc  the  following.  For  the  four-coordinated  tetrahedron, 

4 't*  ,  1  =  10  (3/8) 3/2  ill  [1  -  (3/8)3/2o  /R3]  (2. 

cm  rj  s 
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For  the  s ix -coord inated  octahedron, 
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(2.23) 


(2.24) 


(2.25) 


c.  The  total  potential  energy 


The  total  potential  energy  function  is  obtained  by  adding 
a  repulsion  to  the  sum  of  and  c$...  Thus,  the  total  potential 

function  is  written  as 


CU  (R) 


(2.26) 


For  some  models,  the  repulsion  is  combined  with  a  van  der  Waals 
attraction,  as  is  the  case  for  the  Lennard -Jones  potential.  When  it 
is  clear  that  a  van  der  Waals  attraction  between  the  ion  and  solvent 
is  implicit  in  4>r,  then  specific  account  of  the  interaction  is 
omitted  from  cqn  (2.23)  and  (2.26),  (and  hence  in  (2.26)),  viz.  , 


Rfi 


2c 


^16- 


is  omitted. 

Mechanical  equilibrium  is  specified  in  terms  of  Ro,  the 
equilibrium  value  of  the  radius  of  solvation,  by 

U* (Ro)  =  o.  (2.27) 


Thus,  we  can  write 
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or 
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and 


"  d)  * 


R1?*1 


m  +  n(R0/ p) 


(2.31) 


for  the  Woodcock  generalized  exponential  repulsion. 

Eqn  (2.28)  yields  an  equation  in  the  adjustable  parameters. 

The  complete  expression  for  the  force  constant  yields  another  equation. 
Thus,  with  the  assumption  of  the  value  of  the  solvation  radius 
at  equilibrium  and  the  value  of  the  force  constant  for  the  far 
in frared - act i ve  vibrations,  it  is  possible  to  determine  the 
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values  of  the  adjustable  parameters  for  a  particular  ion.  'this 
process  is  carried  out  in  the  next  section. 

3.  Analyses  with  specific  functional  forms 

The  analyses  to  he  considered  conviently  divide  into  those 
which  depend  upon  exponential  forms  for  the  repulsion,  and  those 
which  depend  upon  inverse  powers. 


The  simple  Horn  exponential  repulsion  is  contained  as  a 
limit  in  the  Woodcock  generalized  exponential  repulsion.  Thus, 
to  begin,  we  consider  the  Woodcock  form. 

Following  the  prescription  laid  down  in  the  last  section, 
we  proceed  here  to  determine  the  values  of  the  adjustable  parameters. 
From  the  equilibrium  condition,  eqn  (2.28),  we  can  write  the 
fol low ing 

(  7  j 

,  ,  a  e  *  a .  u 

B(m  +  nx)  =  R™  2—  +  2-5 —  +  1  2— -  + 

k  Ro  Ro  Ro 

=  Pc  f  s .  n 

in  which  we  have  defined  x  as 

X  =  ( R  0  /  P ) n  .  (3.2) 


The  expression  for  the  force  constant  can  he  expressed  as 
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B[m(m-1)  -  n(n-2m+l)x  +  n2x2] 


3R 
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l——— 
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+  20 


(3.3) 


These  two  equations  can  be  solved  for  B  and  x.  The  elimination  of 
B  yields 


(n-2m+l)P  +  Q  , 

x  =  - - - -  ±  — 


2nP 


2nP 


[n2+2n-4mn-4m+5] P2  +  2(n+l)PcQc 


1/2 


(3.4) 


The  value  of  p  is 


P 


Ro/x 


1/n 


(3.5) 


Finally,  the  value  of  B  is  found  to  be 

B  =  P  / (m  +  nx) .  (3.6) 

The  Born  repulsion  follows  from  the  Woodcock  form  with  m  =  0 
and  n  =  1.  Values  of  B  and  p  for  the  Born  potential  are  given  in 
Table  2. 

4a 

Woodcock^  in  studying  the  alkali  halides,  assumed  a  value 
of  m=* 4 .  The  value  of  n  ranged  from  1  to  6.  We  assume  values  of 
m  =  4  and  n  =  6.  The  results,  values  of  B  and  p  for  the  solvated 

alkali  cations,  are  given  in  Table  3. 


-19- 

l [ ,  Inverse  powers 

lip  to  this  point  we  have  used  forms  for  the  potential  which 

depend  upon  two  parameters  within  the  repulsion.  It  is  useful  and 

instructive  to  shift  attention  from  these  exponential  forms  in  order 

- 12 

to  consider  another  common  iorm  of  interaction,  the  R  repulsion. 
When  this  form  of  repulsion  is  combined  with  an  R  ^  attraction  as 

<K,  =  A/R12  -  B/R6  (3.7) 

M 

1  1 

we  have  an  example  of  the  Mie  potential.  "  On  the  other  hand,  these 
two  terms  can  be  combined  as 


*LlJ  =  4Df  (o/R)12  -  (o/R)6] 


l  3  .  S  ) 


which  is  a  Lennard -.Jones  potential J ^  The  quantity  o  is  a  distance 
for  which  the  potential  vanishes.  F:or  our  purposes,  it  is  more 
convenient  to  use  (3.8)  in  the  form 

*LJ  =  e(p/R)61(p/R)6  -  2]  (3 

L  =  -D. 

When  eqn  (3.7)  or  (3.9)  is  combined  with  other  terms  to  give  a 
potential  of  the  form 


U  (R)  " 


a  e 
„  s 

-  c — -  -  c - 

R2  2R4 


^dd  +  ^ (M  or  LJ) 


I3.1()i 


tlic  quantities  A  and  B  or  c  and  p  arc  parameters  which 
arc  to  he  determined  with  the  use  of  the  experimental  data. 
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To  begin,  we  examine  the  use  of  (3.10)  with  the  Lennard -Jones 

component.  This  form  of  potential  is  analogous  to  the  Stockmayer 

form  of  hybrid  potential  which  is  used  in  the  analysis  of  polar 
14 

molecules.  As  we  illustrate  momentarily,  the  use  of  this  form 
of  potential  with  the  values  of  the  quantities  R0,  a  ,  and  p  in 
Table  1  leads  to  a  mathematical  pathology;  an  intrinsically  positive 
definite  quantity  is  determined  to  be  negative. 

The  complete  force  constant  for  the  potential  with  a  Lennard- 
Jones  component  is 

ct  c  ^ 

k  =  4c-S-  (p/R0)6[11(p/Ro)6  -  5]  -  2c— —  .  (3.11) 

Ro  Rjj 

Define  the  quantity  x  as 

x  =  (p/Ro)  6  •  (.3.12) 

Then,  with  the  use  of  U'(R0)  =  0,  we  find 

ex(x  -  1)  =  A 

(3.13) 

ex ( 1 1 x  -  5)  =  B 


with 
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and 
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The  pair  of  equations  13.13)  are  solved  easily  for  x  and  c  : 
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Here,  the  pathology  appears.  For  lithium,  or  any  of  the  cations 
of  group  I,  we  find  x  <  0.  By  definition,  however,  cqn  (3.12),  x 
must  be  greater  than  zero.  Hence,  the  contradiction  suggests  either 
[1J  the  use  of  the  hennard- Jones  component  is  inappropriate  or  (.2) 
the  values  of  R0  used  are  inappropriate.  We  return  to  consider  this 
question  briefly  in  the  next  section. 

flic  Mic  potential  is  a  decoupled  form  of  Lennard- Jones  potential. 
The  adjustable  parameters  A  and  B  can  be  determined  easily  as  follows. 
The  complete  force  constant  is 
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b  =  B/Ro 

Then,  wo  find  the  equations 
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Thus , 


a  =  pjO  -  5X/5) 
b  =  g-(Y  -  1 IX/ 3)  . 
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13. 


The  results  of  the  application  of  these  equations  to  the  alkali 
metal  cations  are  given  in  Table  4.  It  is  remarkable  to  note  that 
for  both  lithium  and  sodium,  the  value  of  B  found  is  negative .  Thus, 
for  both  of  these  ions  we  find  an  R  ^  l  uleion  in  addition  to  the 
R  ^  repulsion  when  c  =  4. 


4.  Discussion 


I;or  some  time  now,  the  problem  of  the  determination  of  the 
energies  of  solvation  has  commanded  much  attention.*^  These 

efforts  amply  show  that  electrostatic  contributions  predominate. 

Some  calculations  have  been  carried  out  which  incorporate  an 
accounting  for  repulsive  interactions.^’^'  As  one  expects,  the 
repulsion  is  only  a  small  contribution  to  the  total  energy  of  solvation. 
Thus,  one  would  anticipate  a  lack  of  enthusiasm  for  the  pursuit  of 
the  most  accurate  form  of  the  repulsion,  if  the  goal  were  only  to 
determine  the  energies  of  solvation. 

As  we  have  shown,  certain  phenomena,  principally  the  far  infra¬ 


red  spectra,  depend  sensitively  upon  the  repulsion.  In  this  paper 
we  also  have  shown  that  these  vibrations  depend  substantially  upon 


van  dor  Waals  and  induced  dipolar  attractions  which  operate  between 
the  ion  and  the  molecules  of  solvent.  The  energies  of  solvation  are 
thermodynamic  quantities  which  are  characteristic  of  a  state  of 
equilibrium.  On  the  other  hand,  the  vibrational  spectra  of  the 
solvated  ions  reflect  the  importance  of  terms  in  the  potential  energy 
function  which  intluence  t i me  -  dependent  phenomena.  Thus,  although 


the  repulsion  and  the  induced  dipolar  and  van  der  Waals  interactions 

i 

i 

may  play  smaller  roles  in  the  determination  of  the  energy  of  solvation,  j 

their  role  can  become  much  enhanced  for  rate  and  transport  processes.  j 

It  is  possible,  of  course,  to  generate  many  different  forms 
for  the  repulsion  in  the  total  potential.  In  the  absence  of  any 
other  measurable  quantity,  there  is  little  upon  which  to  base  a 
choice  of  one  form  over  another.  liven  though  the  repulsion  and 
van  der  Waals  attraction  contribute  but  a  fraction  to  the  total 
energy  of  solvation,  the  magnitude  of  *iie  contribution  for  one 
form  as  compared  to  another  can  be  diagnostic. 

In  the  previous  section  we  considered  three  model  potential 
energy  functions  in  which  only  the  repulsion  contained  adjustable 
parameters.  In  contrast,  for  the  Mic  potential,  we  determined 
coefficients  for  both  the  repulsion  and  the  van  der  Waals  terms. 

Thus,  it  is  necessary  to  compare  only  functions  which  have  the 
same  overall  form. 

Table  5  lists  the  values  of  the  repulsions  and  the  Mie  potential 
evaluated  with  the  equilibrium  radius  of  solvation.  The  Born  and 
Woodcock  terms  all  can  be  compared  one  to  another. 

It  is  clear  that  there  is  little  to  distinguish  the  Born  from  the 
Woodcock  4/6  potentials  for  the  ions  lithium,  sodium  and  potassium. 


i 
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Thc  Born  repulsion  decreases  monoton ical ly .  This  is  true  also  for 

the  Woodcock  potential.  However,  for  the  Woodcock  4/6  potential 

we  see  that  changing  from  4-  to  6-coordination  at  Kb  shows  a  dramatic 

change  in  the  repulsion.  This  change  is  not  true  of  the  Born 

1  7 

potential.  As  Mui rhead- Gould  and  Laidler  have  shown,  there  is 
a  distinct  difference  in  the  trend  in  the  energies  of  solvation 
upon  changing  from  4-  to  6-coordination.  This  trend  is  matched  by 
the  experimental  points.  However,  it  is  difficult  to  say  at  this 
point  whether  the  trend  depends  significantly  upon  changes  in  the 
repulsion  or  not. 

The  Mie  potential  deserves  special  attention.  First  of  all, 
the  parameters  A  and  B  apply  respectively  to  a  formal  repulsion  and 
an  attraction  which  operates  between  the  ion  and  the  molecules  of 
solvent.  In  contrast,  in  the  other  models  of  the  potential  function, 
all  the  adjustable  parameters  apply  to  one  part  of  the  complete 
function,  the  repulsion.  The  use  of  the  Mie  potential,  therefore, 
suggests  that  the  representation  is  more  physically  ’’truthful”  than 
is  the  case  for  the  other  more  restricted  models.  There  is  no 
compelling  reason  to  expect  that  the  experimentally  measured  (or 
more  accurately,  estimated^)  polarizabilities  for  the  ions  should 
apply  to  the  van  der  Waals  interaction  between  the  ion  and  the 
solvent.  The  results  with  the  use  of  the  Mie  potential  seem  to 
indicate  that  the  experimental  pol ar i zab i 1 i t i es  of  the  ions--as 
used  for  example  with  the  exponential  forms  in  section  3,  part  I-- 
may  indeed  not  apply. 
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The  values  for  the  contribution  of  the  Mie  potential  to  the 
total  energy  of  solvation  are  listed  in  Table  5.  These  values 
represent  the  combination  of  the  repulsion  and  the  attraction. 

It  is  clear  that  the  total  effect  of  both  terms  is  a  net  repulsion 
at  the  equilibrium  radius  of  solvation.  The  fact  that  there  is  a 
trend  downward  from  lithium  to  sodium  and  back  upward  from  sodium 
to  caesium  is  accounted  for  bv  referring  to  the  changing  role  of 
the  R  ^  term  as  shown  in  Table  4.  The  ionic  polarizabilities  are 
given  in  Table  1.  For  lithium  and  sodium,  these  quantities  are 
small.  Thus,  any  induced  dipolar,  van  der  Waals  attraction  between 
the  ion  and  a  molecule  of  solvent  will  be  extremely  small.  For 
example,  for  lithium  we  calculate  this  contribution  to  be  77  J  mole 
It  would  come  as  no  surprise  therefore  to  find  small,  but  positive, 
values  of  B  for  lithium  and  sodium.  It  is  somewhat  surprising  to 
find  in  fact  that  the  polarizability  has  vanished  (in  terms  of  this 
model)  and  that  the  van  der  Waals  attraction  is  replaced  by  a 
repulsion.  Quantum  mechanical  calculations  to  determine  the 
nature  of  the  interaction  between  the  lithium  cation  and  the 
formaldehyde  molecule  indicated  that  the  lithium ^cati on  does  not 
participate  in  any  covalent  delocalization  of  electrons  in  either 
the  Is  or  2s  shells.  That  work*^  concluded  that  lithium  behaves 
essentially  as  a  hard,  charged  sphere.  In  light  of  this  finding, 
the  role  reversal  of  the  R  ^  term  looses  its  surprisal  value. 
Lithium,  and  to  a  much  lesser  degree  sodium,  is  simply  a  very 
classical,  charged  particle. 

It  is  of  interest  to  sec  if  there  is  any  observable  relation- 
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ship  between  the  values  for  the  B-factors  which  are  determined  hv 
matching  the  potential  to  the  experimental  quantities  and  the  values 
of  the  polarizabilities  (listed  in  Table  1).  IV e  have  seen  above 
that  for  lithium  and  sodium  there  can  be  no  correspondence.  The 
fact  that  a  negative  value  of  B  is  found  both  for  lithium  and 
sodium  cannot  support  any  interpretation  which  involves  the 
polarizability  of  the  ion  (viz.,  a  van  der  Waals  type  of  inter¬ 
action).  For  the  ions  of  potassium,  rubidium,  and  caesium  we  find 
positive  values  for  the  R-quant it ies ,  as  given  in  Table  1  . 

At  first  glance,  it  does  appear  that  there  is  some  correspondence 
between  the  B-values  and  the  experimental  values  of  the  polariza¬ 
bilities.  We  can  assume  that  B  is  related  to  the  polarizability 

4 

of  the  ion  and  the  dipole  moment  of  the  solvent  by 

B  =  fot^p 2  .  (4.1) 

The  quantity  f  is  a  constant  of  proport i  onal  i  tv .  Table  (>  lists 

the  values  of  feu  with  f  =  1.  All  the  values  listed  are  very 

large,  much  larger  than  the  experimental  polarizabilities.  If  we 

assume  that  these  values  are  large  because  f  should  not  be  unity,  then 

a  constant  value  of  f  should  emerge  when  we  examine  the  ratio 

a.r  i  ./a.r  ..  These  ratios  arc  also  included  in  Table  b.  As  is 

i  (calc)  l  (exp) 

readily  seen,  there  is  no  cons i stancy .  Lastly,  we  list  the  values 
of  ai(n)/ai(n  1)  **or  calculated  and  the  experimental 

quantities.  There  is  no  similarity  here  either. 

We  conclude  that  for  the  Mic  potential  at  least,  there  is  the 
indication  that  the  trends  in  the  polarizabilities  of  the  ions 
are  roughly  followed.  However,  because  the  effective  polarizabilities 


are  much  larger  than  the  experimental  quantities,  there  appear  to 
lie  operating  additional  quantal  effects  [which  would  involve 
electronic  overlap]  which  are  not  adequately  accounted  for  in  terms 
of  the  simplicity  of  the  potential  energy  functions  used. 

The  results  obtained  for  the  Mie  potential  suggest  that  the 
pathology  observed  with  the  use  of  the  Lennard- Jones  component  is 
intrinsically  associated  with  the  form  and  not  with  the  values  of 
the  parameters  used.  The  Lennard -Jones  potential,  when  used  by 
itself,  adequately  represents  the  balance  between  attraction  and 
repulsion  with  two  parameters,  an  energy  and  a  distance.  When  used 
in  conjunction  with  a  strong  interaction  such  as  an  electrostatic 
ionic-dipolar  interaction,  it  may  not  be  possible  for  the  parameters 
rationally  to  adjust  in  a  physically  reasonable  sense.  Thus,  we 
suspect  that  it  is  more  appropriate  to  use  forms  for  the  potential 
which  consist  of  independently  variable  repulsions  and  attractions. 

We  cannot  claim  to  have  unearthed  the  accurate  form  of  the 
potential  energy  function  which  operates  to  describe  ionic  solvation 
However,  we  claim  to  have  accented  some  factors  which  have 
received  less  attention  than  they  warrant,  namely,  the  roles  played 
by  the  repulsion  and  dispersion  forces.  The  question  must  therefore 
be,  what  is  the  accurate  form  of  the  potential.  A  considerable 
body  of  work  in  related  areas^’^  suggests  that  accurate  potential 
energy  functions  can  be  constructed  in  terms  of  individual  atomic 
and  electronic  (bonding  and  lone  pair)  contributions  from  all  the 
atoms  in  the  molecular  aggregate.  Schcraga's  efforts, in  particul 
have  been  notably  successful  in  a  programme  to  parameterize  a 
large  number  of  interactions.  A  feature  of  Schcraga's  work  which 
seems  also  to  be  applicable  to  the  solvated  ion  is  the  use  of  a 
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three-paramet er  potential  for  each  "electron”.  The  electrons  in 

Scheraga’s  model  are  negative  point  charges  which  are  not 

necessarily  located  at  the  same  positions  as  the  actual  quantal 

J () 

electronic  charge  distributions.  The  potential  Scheraga  ‘  uses  for 
the  interaction  between  electrons  is  the  Buckingham  exp-6  form. 
Pispersive  interactions  are  associated  only  with  the  electrons. 

The  parameters  are  then  fixed  by  means  of  least  squares  fits  to  all 
available  experimental  data  and  to  some  quantum  mechanical  calcula¬ 
tions  in  regions  where  strong  repulsions  operate. 

- 1  2 

Unless  one  is  content  to  use  the  R  or  similar  form  ot 
repulsion  in  the  potential  energy,  then  the  use  of  an  exponential 
form  for  the  repulsion  requires  two  adjustable  parameters. 

As  we  have  seen  with  the  Mie  potential,  it  is  useful  to  he  able 
to  determine  the  coefficient  of  the  van  dcr  Waals  term  by  fitting 
to  experimental  data  as  well.  Thus,  to  be  able  to  fit  both  the 
repulsions  and  the  dispersive  forces  requires  a  minimum  of  three 
adjustable  parameters . 

The  use  of  the  exponential  form  for  the  repulsion  has  a  long 
tradition.  And,  its  use  seems  to  be  more  closely  related  to  actual 
quantum  mechanically  calculated  repulsions  than  is  the  case  lor  the 
R'1  repulsion. Thus,  wc  suspect  that  the  accurate  lorm  ol  the 
potential  energy  function  for  a  solvated  ionic  system  which  will 
emerge  will  be  one  which  faithfully  reproduces  the  vibrational 
spectra  of  the  ion  in  the  cage  of  solvent,  and  the  cage  itsell. 
Moreover,  the  function  ought  to  be  able  to  predict  vibrational 
contributions  to  the  activation  of  rate  and  transport  processes. 

It  seems  at  this  time  that  a  hybrid  form  which  uses  the  exponential 
repulsion  together  with  the  van  dcr  Waals  attraction  for  part  of 
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the  potential  will  be  successful  . 
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Table  1:  Experimental  data  for  use  in  determining  potential  energy 

functions 


Solvent:  dimethyl  sulphoxide:  n  =  1.4783  (Na  D-line),  y  =  4.3  debye 


Ro(A) 


a  =  8A; 


v (cm 


k (dyn  cm  ^ )  ^ 


(A  3 ) 


3>(d) 


4  -coord 

6-coord 

.  + 

Li 

3.24 

429  +  2 

74307 

74894 

0.08 

Na+ 

3.56 

198  +  3 

49497 

50652 

0.20 

K+ 

3.79 

153+3 

47837 

49672 

0.87 

Rb  + 

3.91 

123  +  3 

60121 

64785 

1.43 

+ 

Cs 

4.12 

110+4 

66506 

73888 

2 . 50 

(a)  Refractive  index  and  dipole  moment:  Merck  Index  (Merck  Publishing 

Co.,  Raway  N.J.,  1976);  polarizability  from  refractive  index 
and  eqn  (2.6) . 

(b)  From  Maxey  and  Popov,  see  ref.  2 

2 

(c)  Calculated  from  Maxey  and  Popov's  data  for  4-  and  6-coordination, 

see  ref.  1. 

(d)  From  Smyth,  ref.  6. 


Table  2:  Model  potential  with  Born  exponential  repulsion 


Ion 

B*1012erg 

pxlO^cm 

4 -coord 

6-coord 

4  -coord 

6  -  coord 

Li  + 

0.51 

0.59 

0.26 

0.32 

Na+ 

0.38 

0.43 

0.27 

0.34 

k+ 

0.27 

0.30 

0.24 

0 . 50 

Rb+ 

0.18 

0.20 

0.18 

0 .22 

+ 

Ls 

0.12 

0.13 

0.15 

0.17 

Table  5: 

Model  potential  with 

Woodcock  4/6 

exponential  repulsion 

Ion 

44 

BxlO^erg- 

■  cm1* 

pxi O^cm 

4 -coord 

6-coord 

4  -  coord 

6  -  coo  rd 

Li  + 

59.5  7 

176.53 

3.02 

3.38 

Na  + 

63.18 

186,91 

3 . 29 

3.6  7 

K+ 

59.38 

172.41 

3.24 

3.7  7 

Rb  + 

49.16 

142.05 

3.40 

3.7  2 

Cs  + 

42.62 

122.49 

3.46 

3.7  7 

i 


Table  4:  Model  potential  with  Mie  components 


Ion 

Axio1 

:02erg  cm12 

Bxio 

58erg  cm6 

4-coord 

6-coord 

4 -coord 

6-coord 

.  + 

LI 

0.68 

0.49 

-0.75 

-0.56 

Na+ 

0.69 

1.23 

-0.01 

0 . 56 

K+ 

4.43 

3.20 

9.57 

6.99 

Rb+ 

9.71 

7.22 

31.88 

24.22 

Cs  + 

23.96 

18.20 

71 .23 

55.25 

energy 

Table  5:  Differences  in  4  as  a  contribution  to  the  free/of  cationic 

rep 

solvation 


Ion 

Born 

kJ  mole  * 

Woodcock 

Mie 

Li  + 

30.71 

32.55 

30.88 

Na+ 

22.89 

23.68 

10.04 

K+ 

16.28 

17.91 

29.04 

Rb+ 

12.06 

36.61 

31.21 

Cs  + 

7.82 

25.61 

41.09 

i 


Table  6:  Effective  polarizabilities  from  B  =  oup2 


Ion 

ai  (calc) ^  ^ 

aiCn)/“i(n-l)(a) 

acalc^aexp 

K+ 

J. 

52 

2.52(1.64) 

60 

Rb 

4. 

131 

2.28(1.75) 

92 

Cs 

299 

120 

(a)  Parenthetic  values  are  listed  for  the  experimental  quantities. 

*  *  •  1  "tu. tv  ol  t  o  ion  xn  the  sonso  of  increasing 


atomic  number. 
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